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Abstract: A method is described by means of which, the l ami nar friction 

layer at a wall with arbitrary pressure distribution may be 
calculated from Prandtl's boundary— layer equations. 


Outline : 


I. INTRODUCTION 

II. TRANSFORMATION OF THE BOUNDARY-LAYER EQUATIONS INTO THE 
HEAT-CONDUCTION EQUATION 

III. A SECOND TRANSFORMATION (NEIGHBORHOOD OF THE STAGNATION 
POINT) 

IV. TRANSFORMATION OF THE BOUNDARY-LAYER EQUATIONS INTO 
THE BEAT-CONDUCTION EQUATION ACCORDING TO RRANDTL 
‘AND MISES 

Y. THE NEW CALCULATION METHOD 
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. Symbols 

x,y coordinates parallel and vertical to the wall 

u,v velocity components parallel and vertical to the wall 

U(x) prescribed velocity distribution 

Y stream function 

Uco reference velocity (free— stream velocity) 

L reference length (profile chord or half the profile circum- 

ference, respectively) 

UJL 

Re = -ap Reynolds number 

*"Die allgemeine Losung der Prandtlschen Grenzschichtgleicbungen." 
Lilienthai-Gesellschaft fur Luftfahrtforschung Bericht Ihl, Oct. 19^1, 
PP. 3-7. 
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I. IKTEODUGTIOn 


The solutions of* Prandtl's boundary— layer equations (reference 1) 
known so far (compare Howarth's report (reference 2) and E. Schroder’s 
lecture on the Gottingen boundary— layer meeting) represent special 
solutions inasmuch as they are valid only for special pressure 
distributions and special initial values. For general application, the. 
problem requires the solution of how a velocity profile prescribed at a 
point on a solid wall develops along the wall under the influence of a 
given pressure distribution. 

A transformation of the boundary— layer equations, so far probably 

unknown,- 1 - will be given below with the aid of which the problem may be 
traced back to one already known: the solution of the heat-conduction 

equation. In order to enable, also, calculation of the „ stagnation point 
profile, a second transformation is used which by a wholly analogous 
process also leads to the solution of a linear differential equation. 

The connection between the new transformation and the transformation into 
the heat-conduction equation previously indicated by Prandtl (refer- 
ence 3 ) and Mises (reference k) is shown. 

By means of the resulting calculation method, the la m inar stagnation 
point profile and the boundary layer on the circular cylinder may be 
calculated. 


II.- TRANSFORMATION OF THE BOUNDARY-LAYER EQUATIONS INTO 
THE HEAT— CONDUCTION EQUATION 


In calculating the laminar friction layer on a wing profile for high 
Reynolds numbers, one may presuppose the radius of curvature of the wall 
as large in proportion to the dimensions of the. boundary layer so that 
the wall curvature may be neglected. It enters Into the “calculation 
only indirectly because the influence of the pressure distribution along 
the wall, which Is impressed on the boundary layer from the outside and 
which Is known from measurements or calculation, is considered. 

“'"After conclusion of this report the author's attention was called 
to a report by Piercy and Preston (reference 12) in which a similar 
transformation of the boundary— layer equations is performed for the case 
of a constant pressure distribution. 
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If u and v denote the components of the -velocity in the 
friction layer in x— and y— direction, respectively, that is, parallel 
or vertical to the wall y = 0, V , the kinematic viscosity of the air, 
and TJ(x) , the velocity distribution outside of the friction layer to 
be calculated by means of Bernoulli’s equation, u and v must satisfy 
the force equation 


8u 8u _ dU , 8^u 

u + v — = U — + v — — 

8x 5y dx £y d 


( 1 ) 


and the continuity equation 


hi + h = 0 (a) 

8x 8y 


The latter is to be fulfilled by introduction of the stream 
function “"Kxjy) 


.St, v—St ■ (2a) 


One may put 


*(x,0) = 0 (3a) 


As boundary conditions one has at the vail the condition of no slip 


u(x,0) = 0, v(x,0) = 0 


(3b) 


and for large distance from the wall the end condition 


u(x, 00 ) = U(x) 


(3c) 
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In addition, there has to he a suitable initial condition; u(0,y) must 
he known. If the friction layer begins at a stagnation point (tj(0) = CD , 
the most' important case for the wing, one has to put 


u(0,y) = 0 


w 


It is shown that a linear differential equation may he obtained 
from the quadratic equation (l), by first solving (ll asymptotically 
for large y, thus linearizing with respect to small values of U — u. 
If L denotes a fixed reference length, U M a constant velocity, and 

Re = - 2 ^ f a Reynolds number, one may write (l) also In the form 


Be 


-l/n\~ 1 T 2 d 2 U(U — u) 

■O L 


- L 


L. u ( p - u ) 

Sx Ua, 2 


*\ 


+ 


zS- l — - -B.) 

U Sy 



= L 



u - ^ dy - 2^1 H— -a 

3k U # U Uoo 


IT — u ci_ TJ — u _ ytT* ^ U — u\ 
U„ \ ax U M U 3y U„ J 




(5) 


For large y the terms of the right side are small compared to those of 
the left so that one obtains, for large values of y, a linear equation 

for the function The latter may be transformed Into the heat— 

conduction equation by passing to new coordinates 4, tj. These are 
defined by 


4 = 


TJ(x) dx 
ILL 




tj = JRe 

‘ LDoo V 




(6) 
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thus 


iL = (3L\ 1 ( l s ytt 1 L 8 \ 

\TjJ \ 8x U ffl 3y/ 
Stj V u »/ 


If one finally puts 



and 

5dUl = s(l) (8) 

u. 2 


there results 


= u(u - *) 

^ TJ m 


(9) 


and one obtains from ( 5 ) the differential equation 


a% _ a g sa m as|~a 2 m , 

a ^ 3 si a,' an 2 asW &n as an®-' 
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or written differently 


ih. _ _2fL = u_ fat & _ i£ _a?L 

34 3tj a(4) _34 3 tj 2 3tj 34 3tj 



( 10 ) 


The boundary conditions now read 


£(4,0) = 0, £ (4,0) = a( 4), 


^( 4 ,°°) = 0 


( 11 ) 


and the initial condition for the case of the stagnation point is 


^(OjTj) = 0 


( 12 ) 


The transformation determinant has the value 


proportional _ 
xy so that 
x = 0 

reversible _ 
condition 
disturbance . 

For calculation of the stagnation point profile, however, one must turn 
to other coordinates (compare section III). 

For a general solution of equation (10), one first calculates a 

first approximation of — by putting the left side of (10 ) equal to 

zero and satisfying the boundary and initial conditions (U), (12). 

After introducing this first approximation into the right _ side one 
determines a second approximation from (10) with known right side. 

In this manner, the solution can, step by step, be improved by iteration. 
The essential fact is that the first approximation satisfies both boundary 
conditions, and completely solves equation (10) above all for large t), 
where the solutions contain an essential singularity. 


^X^y ~ ^ y*lx 


IJ^x) y^Re 




thus becomes zero at the stagnation point. 

r2 


If others U _is 

to x, | becomes proportional x^ and tj proportional .to 
the straight line | = 0 corresponds to the straight line 
although the correspondence is not a unique, point by poljxt 
one. Since thereby the boundary conditions and the initial 
are not infringed upon, this singularity does not cause any 
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The iteration method may also he ■written as a series development. 
If denotes the difference of the and (E — l)4h approxi- 

mation for & one has 




$ 

K=1 


•K> 



(13) 


with £ E satisfying the equation (E = 1, 2, . . . ; R]_ = 0) 


St, 2 81 


%( ^Tj) 



with the boundary conditions 


^l( ^0) = a( I) 

^E(M) = 0(E > 2) > 


= 0(E^ 1) j 


and at the stagnation point the initial condition 

£ K (°,Tl) = 0(E > 1) 


(14) 


(15) 


(16) 
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This is understood at once if equation ( Ik ) is summed over K and the 
resulting double sum of the right side rearranged. Therewith the 
solution of the boundary—! layer equation is traced back completely to a 
known problem, the treatment of the inhomogeneous heat-conduction 
equation. 

The shearing stress distribution t in the friction layer 'is 
obtained from 


pu 


s/rs -s-i; 

E=1 Stj 


the displacement thickness 


S* 


-TM) 




from 


g u)5«(x) . . is rk d „ 

v a( § ) K>l(/0 a 


and the momentum thickness 


= r.±(- L _ ^ 

Jo u ) 




p{x)a«(x) _ u(x|<(x1 _ ^ J* 


i-(i) 

Stj \a / 


dt! 


from 
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IH. A SECOND TRANSFORMATION (NEIGHBORHOOD OF THE STAGNATION POINT) 


Thfe singularity of the transformation (6) at the stagnation point 
(U(x) = xU* (0) + . . .) can he avoided, if a further transformation is 
applied to equation (10). 

If one puts 


S - 1 cap 

1 = 

C - «? X 




> 


c = 



J 


(17) 


one obtains altogether a transformation which is regular also for x = 0. 
Moreover the new variables 
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thuS 


§-(!-») (19) 

^ ^ U'(O)^ U cbc 


are more suitable for representation of the laminar profiles than the 
former ones because they are better adapted to the course of the boundary 

layer, the thickness of which increases approximately like v/x. If in 
addition, one puts 


A(F) = a(o) = 1 (20) 

2c 2 6 


the boundary— layer equation reads in these variables 


a, 

aj3 


- - a - X J£L , 1 

=> v— -s — — 




A 




df S^ 2 Sn 81 3n 



(£) 2 > 




- .Ail^.i/ar 21 

Vl 6 <sf 2 W 


(21) 


The boundary conditions are 
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As initial condition T(0,q) lias to "be prescribed so that the ordinary 
differential equation originating from (21) for T = 0 is satisfied 
with the boundary conditions (22). 

Since, with the aid of (17), Green’s function of the linear 
differential equation for the rectangular region, which is found by 
putting the left side of (21) equal to zero, also may be guessed, the 
singular transformation (6) could as well be avoided and the entire 
calculation performed with the aid of equation (21) j however, for the 
numerical calculation the us© of equation (10) is more convenient. 


IV. TRAUSFOmiATION OF THE BOUNDARY— LAYER EQUATIONS INTO 
THE HEAT-CONDUCTION EQUATION ACCORDING 
TO ERANDTL AND MISES 


Erandtl (reference. 3) and Mises (reference 4) have shown that one, 
likewise, obtains a form of boundary— layer equations related to the 
heat— conduction equation, if one introduces, instead of the distance 
from the wall, the stream function \)r as new independent variable 


and. puts 



(23) 


z* = 


U 5 - u g _ 

~~u7~ " 



( 24 ) 


One must then solve the equation 
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with the boundary conditions 




£*(i*,o) = £ 0 * = 


u« 


£*(!*,«) = 0 


(26) 


and a suitably prescribed initial distribution £*( 0 ,ti*). Asymptotically 
it is, for large tj*, also transformed into the heat— conduction equation, 
so that a similar iteration method could be used as for 'equation ( 10 ). 

An indication to that end is to be found in Karmah and Millikan (refer- - 
ence 5) ,* however, the disadvantage of this method is that here a 
singular point of the solution lies at the wall as well (the second 


derivative — there becomes, infinitely large), whereas the solution 

n. . . 

of ( 10 ) is regular at the wall and therefore more appropriate for an 
iteration method. 


V. THE HEW CALCULATION' METHOD 


In order to calculate a laminar friction layer, one has, therefore, 
to solve the system (lk) with the initial and boundary conditions (15)', 
(l6). The solutuions can be represented (compare, for instance, Frank— 
Mises, reference 6, page 872) with the aid of Green’s function 


G(i,Tj;x,y) 



( 27 ) 


for the rectangle 


0 S x < 
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as definite integrals in tlie following manner 


£l(M) =~: J .a(x)G y (6^;x,0) dx 


( 28 ) 


and for K ^ 2 


5^(1 , tj ) “ 

2n/?T ,j Q 


if n “ 

/ J dx ( 2 9) 


Unfortunately, the integration in ( 28 ) can "be carried out generally 
only when a( I) is a polynome in I, The integration ip. ( 29 ) is even 
then not generally performable so that numerical methods must he applied 
for the evaluation of the integrals. 

If one puts 


°1 


Jl r X . 

2 \/g — X 


( 30 ) 


t] + y 

cj 2 = — 

2 v/ |- x 


and denotes 



the Gaussion error function, G- may he written in the form 


G(£,TijX,y) = & 


^(cTg) 
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Hence, one obtains from (28) and (29) the relations 


^l(M) = / a ( x ) d.( 4 (cri)) Y =0 


£>*]) “ 


/x=0 

■*(/• 

/x=o \ uy=o 


(-* K (x,y)) d | (K > 2) 


( 31 ) 


Thus one has to plot a against 



and to integrate graphically, or, respectively., ' to plot for fixed 

first against . - - — * - - 


| (®(ai) -®(ff 2 )_) 


and to integrate and to plot the integrals obtained oncemore against 
and to integrate them. The derivatives required in the calculation of 
the Eg values are obtained from (K ^ 2) 
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(/x=0 


( 


7=0 


r 2 2' 

T 1 ~°2 

(-%(x,y)) d[- + e - - 1 1 dx 

2\Zn\/t -xj 


b 2 £- 


/. 


'K 




/x=0 


17=0 


2 2 > 

(-%(x,y)) d(~ ~— 2 -- / / <3^ 

2 Rl -x) 


sr = 5^~ %UjTl) 


Corresponding formulas are valid for the derivatives of 

Since this calculation has to be carried out for various pairs of 
values i, i^, the procedure is rather extensive even though the abscissa 
functions may be calculated once for all. Since for determination of 
a veloc it y profile at a point I only the profiles lying shortly ahead 
need be known more exactly — the profiles lying farther back having little 
influence on the calculation at the point I — information about a few 
profiles at larger distances before and about several profiles shortly 
ahead is sufficient for calculating the profile at a point i. 

tfo general statements may be made so far concerning the convergence 
of the iteration method or an estimation of errors. The numerical calcu- 
lation (compare section Yl) proved that the tg values rapidly decrease 

in the proximity of the stagnation point whereas in the further course 
of the boundary layer, particularly in the neighborhood of the separation 
point, the convergence seems to be less favorable. 


VI. EXAMPLE 


In order to test the iteration method, at first the velocity 
distribution at, the stagnation point was calculated with the aid of the 
ordinary differential equation resulting from (21) for f = 0. In 
figure 1 , the first approximations are compared with the solution 
calculated by Siemens (reference 7 ) and Hartree (reference 8 ). Thus the 
solution is approximated from one side by the iteration method. 
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The laminar friction layer' on a circular cylinder was calculated 
as a further example. The potential theoretical velocity distribution 
(fig, 2) 


l = 26inH 
U_ L 


not the experimental one according to Hiemenz (reference 7) was taken 
as a basis because the function a( |) may then be written as a 
polynomial in I. 


a( 5 ) 


= U = 4«i - *2 8 2 

Uoc 2 


The integral may then be calculated in closed form. Tlie first 


approximation for — = 1 


^1 


is plotted in figure 3 for the points 


jtf; = 0 (stagnation point), Jt ! =1, ret = 2 (maximum velocity), = 2. 
(location of the laminar separation) and «| = 3 against 


6 



7 

x 


fgz r ux 
V \j 2 J'u dx 



(compare section HI); figure h represents the third approximation 


“= 1 ^ 
u = V 


a 1 


>K 


As a rough estimate shows, the third approximation, in general, still 

yields slightly too large values for the error is growing with 

increasing jtS. In the neighborhood of the wall, the various approxi- 
mations yield partly too large, partly too small values; hence, the 



NACA TM 1278 


17 


point of laminar separation could not be accurately determined for the 
present. The calculation is to be continued in order to obtain more 
accurately the profiles for larger | values as veil. 

A fev velocity profiles at the velocity maximum XT’ = 0 are plotted 
for comparison in figure 5- Thus the profile at the flat plate calcu- 
lated by Blasius (reference 9 ) almost agrees with the third approxi- 
mation. Probably the differences decrease still more if further 
approximations are calculated. One can see that the approximation 
method accor ding to Pohlhausen '(reference 10) also yields relatively 
satisfactory values. In contrast, the profile calculated according to 
an approximation method of the author (reference 11 ) deviates somewhat 
more strongly, even though the result for the displacement thickness 5* 
is practically correct . 

A corresponding comparison for several profiles in the proximity 
of the separation point (fig. 6 ) shows a relatively good agreement of 
the third approximation with the approximation method of the author 
(reference 11), whereas the corresponding Pohlhausen profile and, most 
of all, the Hartree profile deviate more strongly. 

On the basis of these examples and of the cases calculated 
according to reference 11 , where in contrast to the one-parameter method 
of Pohlhausen the previous history of a profile was at least partly taken 
into consideration, the statement is permissible that the previous history 
in the region of the pressure drop behind a stagnation point has little 
influence on the form of a boundary— layer profile, so that a one-parameter 
approximation yields usable values. In the region of pressure rise, 
however, it will — in order to obtain more accurate calculation — 
probably always be necessary to consider the entire previous pressure 
distribution, not only the pressure gradient at the respective point. 


VII. SUMMARY 


The solution of the boundary— layer equations for the laminar flow 
along a wall with prescribed pressure distributions is traced back to 
the solution of a system of linear differential equations of the type 
of the heat conduction equation. Due to the necessarily large 
expenditure in time, this doeB not yet represent a method easily 
applicable in practice; however, one now has the possibility of completely 
calculating a few characteristic cases in order to estimate and sometimes 
improve the usefulness of the known approximation methods. 

Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 4.- A few laminar velocity profiles on. the circular cylinder 

(third approximation). 



Figure 5.- Comparison of a few laminar velocity profiles at the pressure 

minimum (U* = 0). 



Figure 6.- Comparison of a few laminar velocity profiles in the region of 

the laminar separation point. 


NACA-tangley - 6-26-50 -1000 





